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Abstract 



^ ' We use Gaudin's Fermi-Bose mapping operator to calculate exact solutions for the Lieb-Liniger 



model in a linear (constant-force) potential (the constructed exact stationary solutions are referred 



^ , to as the Lieb-Liniger- Airy wave functions). The ground-state properties of the gas in the wedgelike 

Q^l trapping potential are calculated in the strongly interacting regime by using Girardeau's Fermi- 

Cd ' Bose mapping and the pseudopotential approach in the 1/c approximation (c denotes the strength 
of the interaction). We point out that quantum dynamics of Lieb-Liniger wave packets in the linear 

Ch I potential can be calculated by employing an A^-dimensional Fourier transform as in the case of free 

i_Zj' expansion. 
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I. INTRODUCTION 



One-dimensional (ID) Bose gases hold great potential for studying the physics of interact- 
ing quantum many-body systems. In recent years, a growing interest in studies of theoretical 
ID models, first introduced by Lieb and Liniger [l| and Girardeau |2|, has largely been in- 
spired by experimental progress in realizing these models with ultracold atomic gases 3|-[7|. 
In experiments, an effectively one-dimensional system is achieved in elongated and trans- 
versely tight atomic wave guides, loaded with ultracold atoms, where transverse excitations 
are strongly suppressed ffl. These atomic gases are well described by the L.eb-Lirriger 
model [1| - a system of identical Bose particles in ID which interact via ^-function interac- 
tions of strength c. In the limit of infinite interaction strength (c — ?■ oo), the Bose particles 
can be described by the Tonks-Girardeau model [2|, describing an "impenetrable" Bose gas. 
This regime occurs when effective interactions are strong, whereas temperatures and linear 
densities are low J8l-ll0| . In pas sing , we note that the Lieb-Liniger model can also be realized 
in the field of quantum optics ll|. An ubiquitous part of the ultracold atomic experiments 
is the gravitational force, which can be nullified in a precise experiment, where ID atomic 
wave guides are horizontal. However, to the best of our knowledge, exact solutions for the 
Lieb-Liniger gas (consisting of an arbitrary number of particles) in the presence of a constant 
external force have not been constructed yet. 

An important breakthrough in the context of interacting Bose gases was achieved in 1963 
when Lieb and Liniger found exact eigenstates for a ID Bose gas with ^-function (repulsive) 
interactions of arbitrary strength c l|]. The underlying structure of the eigenstates reveals 
the Bethe ansatz, and each eigenstate is determined by a set of A^ quasimomenta, where A^ 
denotes the number of particles. If the quasimomenta obey a particular set of transcendental 
equations [l|], the wave functions will obey periodic boundary conditions. This breakthrough 
was followed by studies of the model with attractive interactions 12|. A superposition of 
Lieb-Liniger eigenstates, where again quasimomenta obey a particular set of transcendental 
equations, can be used to construct the ground state in an infinitely deep box [13[. More 



recent studies of the Lieb-Liniger model (e.g., |14|-|21|) are mainly motivated by the growing 
experimental progress 3|-[7| • It should be noted that all of the aforementioned studies focus 
on a system where an external potential is zero in a given region of the x space, and boundary 
conditions are imposed on the wave function at the border (s) of this region (e.g., periodic 



Q, 



semi-infinite line 



are studies of few-body systems 22, 
in a harmonic trap in Ref. 
were presented in 22 1. 



13|, and infini tely deep box 13|). For other external potentials, there 



231]. Exact solutions were presented for A^ = 2 particles 



23| ; solutions for N = 2 and A^ = 3 particles in a linear potential 



One approach in attempting to broaden the set of known exact solutions is by using 
Gaudin's Fermi-Bose mapping operator (call it Or) 24| . which is applicable also to find 



time-dependent Lieb-Liniger wave functions 24|-|28 



In the limit of the Tonks-Girardeau 
gas c — )■ oo, the Fermi-Bose mapping (which was discovered in 1960 |2|) can be utilized in 



any external potential |2| and for time-dependent problems 29J (for a review, see also 30|). 
For the finite coupling Lieb-Liniger gas (c finite), the method of Gaudin has been shown 
to be valid in the absence of any external potential (i.e., on an infinite line 25|]), and has 



been used to study free expansion from localized initial conditions 26 



271]; in this case the 



time-dependent wave function can be calculated via an A^- dimensional Fourier transform. 
Interestingly, such a transform can be also utilized for a Lieb-Liniger gas reflecting from the 
wall [Ml- However, Gaudin's method (at least in its current form) is not applicable to find 
eigenstates of a Lieb-Liniger gas in generic trapping potentials V{x) (such as the harmonic 
oscillator); technically, this arises because the differential operator Oc does not generally 
commute with such potentials. 

Here, we study the Lieb-Liniger model in the constant-force (linear) potential. Exact 
stationary solutions for this system are constructed (we call these wave functions the Lieb- 
Liniger- Airy states) by employing Gaudin's operator Oc- The construction is enabled by the 
fact that this operator commutes with the linear (constant-force) potential. We calculate 
the ground-state properties of the Lieb-Liniger gas in the wedgelike potential [V{x) = ax 
for X > (a > 0), and oo otherwise] in the strongly interacting regime. This is achieved 
in the Tonks-Girardeau regime and below that regime in 1/c approximation by employing 
the pseudopotential approach [32]. Finally, we point out that the time- dependent Lieb- 
Liniger wave packets in the linear potential can be calculated via an A^-dimensional Fourier 
transform. 



II. LIEB-LINIGER-AIRY STATES 



The Lieb-Liniger model describes A^ identical bosons in one spatial dimension which 
interact via (delta- function) contact potential [l|. In this section we consider this system 
placed in a linear external potential. The stationary Schrodinger equation for the many-body 
wave function ipsixi, . . . ,xn) in such a system is 



Eij 



B 



N 



9V 



B 



dxi 



E 2c5( 

l<i<j<N 



N 






:)i'B + ay^^Xii:B, 



(1) 



j=l ' l<i<j<N i=l 

where c > denotes the strength of the interaction, and a > is the constant external 
force. Solutions of Eq. ([T]) for a single particle (A^ = 1) are the Airy functions. For this 
reason, in what follows, we will call the solutions of Eq. ([1]) for A^ > 1 the Lieb-Liniger- Airy 
(LLA) states (we are interested only in those solutions which decay to zero when x — t- oo). 
The constant force in ultracold atomic experiments can arise from the gravity force (e.g., if 
the one-dimensional atomic wave guides are tilted with respect to gravity). 

In what follows, we will demonstrate that LLA states can be constructed via Gaudin's 



Fermi-Bose transformation 2^. Because of the bosonic symmetry of the wave functions, 
one can consider only the fundamental permutation sector of the coordinate space Ri : Xi < 
X2 < ■ . . < xn- Within this sector, the Schrodinger equation ^ reads 

N 



Eij_ 



B 






(2) 



The interaction term is taken into account as a boundary condition (the so called cusp 
condition), which is imposed upon ipB at the borders of Ri (i.e., when two particles touch 



Q): 



1 



d 



dxj+i 



d 
dxj 



^B = 0. 



(3) 



Xj^l—Xj 

Equation ([2]) holds in all other permutation sectors, whereas the interaction cusp (^ can be 
re-expressed on the borders of other sectors as well. To construct the LLA states we utilize 
Gaudin's Fermi-Bose mapping operator 24l |. 

l<i<j<N L \ J « 



Or 



(4) 



which acts upon an antisymmetric (fermionic) wave function ipp- The wave function ipp must 



obey the Schrodinger equation for noninteracting spinless fermions in the hnear potentiah 



N 



N 



E'4jp = -^^-f + a^Xi'4)F. (5) 

The wave function ipp can be written in the form of Slater determinant with Airy functions 
as entries: 



IpF = tt s 



N 



det M[a-^Xj — a ^E^), 



(6) 



^N 



where E = ^^^^ Ei. 

The LLA states [i.e., solutions of the Schrodinger Eq. ([2]), together with the cusp condi- 
tion ([3])], are given by 

i^B,c = KdciJF, (7) 

where A/'c is the normalization constant. It is known that all wave functions of the form ([7]) 
obey the cusp conditions throughout the configuration space 24l-l26| . To show that ipB,c is 
also a solution of Eq. ([2]), it is sufficient to prove that the following commutators are zero: 



J:^^V^xl^, 



and 



J2^^^^0c 



0; this is sufficient because ipp obeys Eq. ([5]). The 



first commutator is trivially satisfied, and therefore we are left to verify that 



/ J Xj, Wc 



0. 



(8) 



As a first step, we restrict ourselves to the case of two particles, N = 2. By using [xj, d/dxi 



-Sj^i, we have 



xi + a;2,sgn(x2 - Xi) + 



I f d d 



X2, 



_d_ 
dxo 



Xi, 



d 
dxi 



0. (9) 



c \dx2 dxi 

Now we generalize this for any number of particles A^. Let us write the differential operator 
as a = ni<i<i<7v^i,i' where 



B 



«j 



sgn X,- -Xi) + - 



d d 



c V dxj dxi 



(10) 
Wi+i ■ ■ ■ Wm, valid for op- 



A general expression, V, UZi ^i = Efii ^i " " " ^i-i ^. ^i 

erators V and Wi, I = 1, . . . , M, enables us to write the required commutator for the case 

of A^ particles: 



^Xk,dc 



2_^^N-l,N 
i<j 



/ J '^ki ^i,j 



■B,: 



(11) 



Now Eq. (j8]) follows immediately because for any Bij we have 'Ylk^kiBij 
Xi + Xj,Bij = 0, as is verified for the N = 2 case. This completes the proof that the 
wave function ipB,c defined in ([7]) is a solution of Eq. ([1]). 

In this section we have found exact closed form solutions of Eq. ([!]). We point out 
that the eigenstates ([7]) with total energy E are degenerate, because the choice of single 
particle energies Ei for which E = J2i=i ^i is not unique. By superposition of degenerate 
eigenstates ([7]), one can construct eigenstates which are of different mathematical form. In 



22j the authors study Eq. ([T]) for N = 2 and A^ = 3 particles. They constructed solutions 
by introducing a new set of coordinates and separating Eq. ([1]). For N = 2 they separate 
the center of mass and relative motion. Their solution for a give n energy can be written as a 



superposition of eigenstates ([7]). For A^ = 3 the procedure in [22] becomes more cumbersome, 
which clearly points out the advantage of using Fermi-Bose transformation for solving Eq. 

©• 

III. THE LIEB-LINIGER GAS IN A WEDGELIKE POTENTIAL: STRONGLY 
INTERACTING LIMIT 

In this section, we consider the Lieb-Liniger gas in the wedgelike potential defined as 

{X if a; > 0: 
(12) 
oo if X < 0. 

For simplicity, we have fixed the value of the constant force to a = 1. Solutions for any 
other value can be obtained by simple rescaling: x — )■ a^^^x and E — )■ a~'^^^E. 

In order to find the ground state in such a potential, one should find solutions of Eqs. 
(|2]) and ([3]) (assuming we work in the fundamental sector Ri), together with the follow- 
ing boundary condition: ipB,c{xi = 0,X2, ■ ■ ■ ,xj\f) = 0. The first idea that may come to 
mind in attempting to find such a ground state is to utilize Eq. ([7]) as an ansatz, since 
it apparently obeys (I2l) and (I3l), and try to adjust the A^ free parameters Ei such that 
'>pB,c{,xi = 0, X2, . . . , xn) = 0. Namely, such a procedure leads to the solutions for the ground 



Q 



unc- 



states of a Lieb-Liniger gas on the ring [![ , where instead of the ansatz (JTj) with Airy 
tions, one utilizes an ansatz with plane waves, iI)b,c = A/'cOc det^j^^ e*'^^^™ (e.g., see 
and instead of Ej, one adjusts the quasimomenta kj (which have to obey Bethe's equations) 



33|), 



to acquire the proper boundary conditions. However, for this wedge hke potential such a 
hne of reasoning fails. Mathematically, this occurs because the first derivative of the Airy 
function is not simply related to the Airy function itself (whereas a derivative of a plane 
wave is proportional to the plane wave itself). 

Nevertheless, we can find solutions in the form ([7]) in the Tonks-Girardeau limit (c — )■ oo), 
and we can utilize some form of 1/c approximation to find deviations from the Tonks- 
Girardeau ground state for large but finite c. The Tonks-Girardeau ground state is con- 
structed by symmetrizing the Slater determinant of A^ lowest single-particle eigenstates |2|: 

1 ^ 

tpTG = Ilk<mSgn{Xm - Xk)^== dct 0i(xj), (13) 

where 

, , , Ai(x — Ei) ., ^. 

The single-particle energies E^ are such that Ai{—Ei) = [i.e., (piiO) = 0], and the eigenstates 
form an orthonormal set: j^ (j)*{x)(f)j{x)dx = 6ij. The ground-state energy is simply Etg = 
Yli=i-^i- ^s an illustration, in Fig. [T]we display the single-particle density for the Tonks- 
Girardeau ground state (dashed blue line) comprising A^ = 10 particles. 

An approximative perturbative approach for calculating the properties of a Lieb-Liniger 



gas in the strongly interacting regime has been suggested by Sen [32|. It can be shown 
that the perturbation around c = oo (the Tonks-Girardeau limit) is correctly described by 
a pseudopotential [32 1 

V;p = -^$^5"(a;.-x,), (15) 

i<j 

that is, the pseudopotential (fTSl) is utilized as a small perturbation around the Tonks- 
Girardeau ground state (unperturbed state) for large c. It gives the correct first-order 
correction to the ground-state energy and wave function when plugged into the standard 
perturbation expressions with 1/c as a small parameter. 

In the 1/c approximation, the ground-state energy of the Lieb-Liniger system is 



EB,c = ETG+(ij- 



'Jtg 



Vpp 



i,^ 



TG 



= Etg--N{N-1). (16) 

Result (TT6|) is obtained by a direct calculation of the expectation value of the pseudopotential 
Vpp for the Tonks-Girardeau ground state. Such matrix elements are readily evaluated by 




FIG. 1: (Color online) The single particle density pB,c{x) (solid black line) of A^ = 10 Lieb-Liniger 
bosons in a wedgelike potential (c = 40, a = 1). Dashed blue line shows the density in the 
Tonks-Girardeau limit. 



using Slater-Condon rules: 



^P^ 



TG 



Vr 



PP 



i,^ 



TG 



V / dxU*{x)(t)i{x)—^[(t)*{y)(pj{y)]y=^ 



'dy' 



0*(x)0j(x)3^[0*(y)(/>i(y)]j;=x 



(17) 



We have verified ( TT6l) numerically (by employing Mathematica) up to A^ = 20 particles, 
and we conjecture that the expression is valid for any number of particles trapped by the 
potential f lT2|) . 



To first order in 1/c, the Lieb-Liniger wave function is given by 32 1 
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Wtg 1 



(18) 



where V'tg" labels an excited Tonks-Girardeau state; this state is obtained from the 
ground state ipTG by replacing the single-particle state 0„,, where n < N, with the single- 



particle state (pm of higher energy, m > N > n. Analogously, iprpQ 



{m,m';n,n') 



labels two 



particle excitation of the TG gas state. The expression for the single-particle density 



Pb,c{x) = N Jdx2 ■ ■ ■ dxN\'ipB,c\'^ can be calculated straightforwardly by employing the wave 
function from Eq. flTSj) . by keeping the terms up to 1/c: 



PbA^) ~ PTG[X) 



Ptg[x) 



N 



n<N,m>N 



,4T^ 



Vr 



PP 



i^TG 



En — E„ 



dxo 



■ dXN ^TGi^TcT^ 
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-,Am;n) 
ypp 



n<N,m>N 



En — En 



■ (f>nix)(f)ni{x) 



Here, the matrix element Vpp 



{m;n) 



Vtg 



E.<, 5"{^^ 






(19) 



ipTG / of the single-particle 



excitation from the level n < N with energy E^, to the level m > N with energy E^, is 
given by 

N 



y{m;n) 
PP 



E 



d"^ d^ 

dx{(pl^{x)(pn{x)j^[<P*{y)(piiy)]y=:r-<Pl^{x)<Pi{x) — [(f)*{y)(f)n{y)]y=x 



(20) 

In Fig. [T]we illustrate the single-particle density Pb,c{x) in 1/c approximation (solid black 
line), which is obtained by using Eq. ( ITQl) for A^ = 10 and c = 40. It should be mentioned 
that the two-particle excitations [second sum in Eq. flTSl) ] do not yield any contribution 
to the first-order single particle density Pb,c{x)i due to the vanishing of the overlap of the 
wave functions in calculation of the density (in the same way as demonstrated for the case of 
bosons confined in an infinitely deep box 32] )• In our calculation of the density ps^c via (fT9l) . 
we have included only a finite number of terms, where the cutoff is chosen to be sufficiently 
large, such that the contribution of the remaining terms is negligible [for the calculation 
illustrated in Fig. [H we kept 150 terms in Eq. ( TT9|) with the highest contribution]. 

IV. EXACT QUANTUM DYNAMICS VIA A FOURIER TRANSFORM 



In this section we discuss the time-dependent solutions of the Lieb-Liniger system in a 
linear potential. Before proceeding, we note that dynamics in the strongly interacting regime 
(i.e., dynamics of a Tonks- Girardeau gas in a linear potential) was studied in Ref. 3J]. Here, 
we assume that the bosons are initially localized by some external trapping potential. At 
time t = 0, this potential is suddenly turned off, and bosons are released to evolve in the 
linear potential. This problem can be related to free expansion of the Lieb-Liniger wave 
packet by simple rescaling of the coordinates. If the wave function ipfree{xi, . . . ,xiy,t) obeys 



the equation, 



.dljjfr 



dt 



N 






(i.e., ipjree desciibes free expansion [26|, |27|), then the wave function 



+ E 2c5( 



Xi 



Xj)^fr 



i)B,c{Xl,...,XN,t) 



-«a<E«=i(a;«+«iV3) 



'^Sree{Xi + Ot^, . . . , Xat + Ot^, t) 



is the solution of the time-dependent problem in the constant-force potential, 
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.di>_ 



B,c 



dt 






dxj 



y^ 2c 5(Xi - Xj)tpB,c + " XI ^i^S.' 
l<i<j<Ar 



(21) 



(22) 



(23) 



The initial conditions coincide (i.e., at t = we have ■ipB,c = '4' free)- Note that the phase 
factor in Eq. (122|) accounts for the momentum per particle at, which is acquired in time 
in the field of constant force a (in units used here, m = 1/2, and therefore the classical 
acceleration is 2a). Transformation (|22l) can be verified by direct substitution in Eq. (|23l) . 
from which it becomes evident that it is valid for any two-particle interaction V{xi — Xj). 
Namely, transformation Xj — )■ Xj + at"^ does not affect the two-particle interaction term 
V{xi — Xj) [in fact, because of this, Eq. (122|) can be deduced from the well-known solution 
for a single-particle wave packet in a linear potential]. 

It is known that freely expanding Lieb-Liniger wave packets can be calculated by solving 
an A^- dimensional Fourier transform 1261: 



^-'-■■■■-"^/— '^ ■-)^-^-''--'^"- 



(24) 



We note that the function G is not the Fourier transform of the wave function ibfree because 

n n 

it depends on the coordinates Xj through the sgn(xj — Xj) terms (see Refs. [26|, |27j for 
details), that is, it differs from one permutation sector in x space to the next. Nevertheless, 
by calculating the integral in Eq. IHM in one sector (say Ri), we obtain ipfree in that sector, 
which is sufficient due to bosonic symmetry. The function G contains all information on 



initial conditions and it can be expressed in terms of the projections o 
function on the Lieb-Liniger free space eigenstates (e.g., see Refs. 



26 



the initial wave 



27(1). By using Eqs. 



21]) and fl22|l we can express ipB,c in terms of an A^-dimensional Fourier transform: 
i/jbA^i^ ■■■^^Njt) = / dki--- dkN 

{ki - at)^ - kf 



N 



X G(/ci,...,/cAr)exp < i^ 



i=l 



{ki - at)xi + 



3a 



(25) 



10 



We would like to note that result (125|) can be obtained by straightforward use of Fermi- 
Bose transformation. The time- dependent wave function ipp which describes the system 
of A^ noninteracting fermions in a linear potential V{x) = ax can be written via its Airy 
transform: 

N 



iPf{xi, ...,XN,t)= dEi--- dENlpF{Ei, ..., ^Ar)e 



-itj:tiE. 



Y[Ai{a 



-2/3, 



aXi-Ei)). (26) 



i=l 



Here, ippiEi, . . . , Ej^i) contains information on initial conditions, 



ME: 



1) 
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TV 



N 



dXi ■ ■ ■ dX]y1pp{Xi, . . . ,X]y 



0)11 Ai 



i(a-2/3(ax,-E,)). (27) 



j=i 



By using the well-known relation between the Airy and Fourier transform ipp 35 1 



MEu 
we find 

i:F{xi,...,XN,t) -- 



^E 



N 



a 



2/3 



N 



rfA;i---rfA;^^^e*S-M-'=«'/=^)/°, 



(28) 



/ dki ■ ■ ■ dkN ipF exp < i \^ 



{ki - at)xi + 



{ki - atf - kf 



3a 



(29) 



The time-dependent solution of the Lieb-Liniger model [i.e. Eq. ( 125|) ]. can now be found 
directly from the expression above by applying the Fermi-Bose transformation operator Oc 
onto Eq. ()29|). 

Our discussion in this section adds upon the previous studies of Lieb-Liniger wave-packet 



dynamics on an infinite line 24N27J]. and in the presence of the hard- wall potential 311; in 



all these cases the motion of an interacting Lieb-Liniger wave packet can be calculated by 
using an A^-dimensional Fourier transform. 

In order to illustrate the connection between (1211) and ( l22i) . we present the following 
numerical example. The system of three Lieb-Liniger bosons are trapped in the ground 
state of an infinitely deep box of length L = tt; at t = 0, the trap is turned off and the 
bosons start to experience the constant force a = 3. The exact initial wave function is 
constructed as a superposition of free space eigenstates [l3|. From this state we can find 
the function G{ki, . . . ,kN) which keeps all information on initial conditions [26|, |27|. By 
numerically calculating the integral in f l25|) . we obtain the time-dependent wave function 
'^B,c(3^i, • • • , a; AT, t) describing the system. Here, we plot two relevant physical quantities, the 



11 



single-particle density, Pb,c{,x, t) = N f dx2 ■ ■ ■ dxN \4'b,c{x, . . . ,XN,t)\ , and the momentum 
distribution n{k,t) (density in k space). 

From Eq. (l22l) it follows that the density in coordinate space will be the same as in 
the case of free expansion {a = 0), with mere translation of the coordinates [pB,c{x,t) = 
Pfree{x + at'^,t)]. The momentum distribution will be equivalent also up to the simple 
transformation k ^ k — at. The density profile and momentum distribution of the wave 
packet are plotted in Figs. [2] and [3|, respectively, for three various interactions strengths c: 
(a) c = 0.25, (b) c = 3, and (c) c = 10. Starting from t = 0, the wave packet evolves to the 
left in X space with the center of mass motion at"^, while at the same time it spreads in width 
independently. For large c, the spread is more pronounced, as can also be conjectured from 
the initial momentum distribution. For very large c, the wave packet will asymptotically 
experience fermionization of the momentum distribution J36|, |37| . 

V. CONCLUSION 

We have studied the Lieb-Liniger model in the constant-force (linear) potential. Exact 
stationary solutions for this system, referred to as the Lieb-Liniger-Airy states, were con- 
structed by employing Gaudin's Fermi-Bose mapping operator Oc- This was enabled by the 
fact that the operator commutes with the linear potential: [Oc, X]„- axj] = 0. We have calcu- 
lated the ground-state properties of the Lieb-Liniger gas, in the strongly interacting regime, 
in the wedgelike potential: V{x) = ax for x > (a > 0), and V{x) = oo for x < 0. This 
was achieved in the Tonks-Girardeau regime and in 1/c approximation by employing the 



pseudopotential approach [32|. Finally, we have pointed out that the time-dependent Lieb- 
Liniger wave packets in the linear potential can be found by employing an A^-dimensional 
Fourier transform. 
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interaction strengths c: (a) c = 0.25, (b) c = 3, and (c) c = 10. Red dotted lines are for t = 0, 
solid black lines are for t = 1, and blue dashed lines are for t = 2. 
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